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ABSTRACT 

The k-filtering technique and wave polarization analysis are applied to Cluster magnetic field data 
to study plasma turbulence at the scale of the ion gyroradius in the fast solar wind. Waves are found 
propagating in directions nearly perpendicular to the background magnetic field at such scales. The 
frequencies of these waves in the solar wind frame are much smaller than the proton gyro-frequency. 
After the wave vector k is determined at each spacecraft frequency fsc, wave polarization property is 
analyzed in the plane perpendicular to k. Magnetic fluctuations have 5B± > (5_B|| (here the || and _L 
refer to the background magnetic field Bq). The wave magnetic field has right-handed polarization 
at propagation angles 0kB < 90° and > 90°. The magnetic field in the plane perpendicular to Bq 
however has no clear sense of a dominant polarization but local rotations. We discuss the merits and 
limitations of linear kinetic Aflven waves (KAWs) and coherent Alfven vortices in the interpretation 
of the data. We suggest that the fast solar wind turbulence may be populated with KAWs, small 
scale current sheets and Alfven vortices at ion kinetic scales. 
Subject headings: solar wind — turbulence — waves 



1. INTRODUCTION 

The solar wind is a natural laboratory to investigate 
plasma turbulence. It is well known that in the inertial 
range, at which the usual magnetohydrodynamic (MHD) 
description is still valid, magnetic turbulence is strongly 
anisotropic: for a given wave number k, magnetic fluc- 
tuation energy is much more concentrated at quasi 
perpendicular propagation {k± ^ fc||) than it is at quasi 

parallel propagation (fc i ^ fcii ) (iShebalin et 1 198 

Goldreich fc Sridhar I 119951: IM atthaeus et al. ' nm- 

2005: Dasso et al. 



Bieber et al. 1 119961 : iHorburv et al. 
20051) . Numerous measurements find the Kolmogorov 



^-5/3 spectrum of magnetic field fiuctuations in the 
inertial range and a steeper spectrum at ion kinetic 
scales (which is often called the dissipation range where 
MHD description breaks down). A spectral break 
point around kpi « 1 (where pi is the ion thermal 
gyroradius) or kdi ~ 1 (where di is the ion iner- 
tial length), which marks the end of the iner- 
tial range, suggests possible initiation of kinetic dissi- 
pation processes at ion scales while turbulent cascade 
continues to operate at the same scales and at smaller 
scales, up to electron gyroradius (Alexandrova et al. 
2009, 2012). It is an open question exactly which 
scale is responsible for the s pectral break (see a rece nt 
discussion on this topic bv iBourouaine et al.l ()2012D '). 
A view to account for the observed spectral steepen- 
ing at high frequencies (ion scales) is to interpret the 
spectral steepen i ng as evidence of kineti c Alfven waves 
(iLeamon et al. I 119981: iBale et al.l 120051: iHowes et al. " 
20081 : iSchekochihin et al. I 120091 : IHowes fc Quataert 



whistler waves (iBiskamp et al. I 119961: ILi et al. 1 120011: 



iSahraoui et al. I l2010al: iSalem et al. I 120121) . or 



IStawicki et aril200lHGarv fc Smith \\20m under the as- 
sumption that although linear waves are unable to pro- 
duce nonlinear cascade, they may still approximately de- 



scribe the nature of turbulence at ion kinetic scales. An 
alternative view is that 2D structures (such as current 
sheets, coherent magnetic vortices) populate the fluctua- 
tions at these scales and have been observed in the iono- 
sphere, magnetosphere and magnetosheath (Chmyrev et 
al. 1988, Volwek et al. 1996, Sundkvist et al. 2005, 
Alexandrova et al. 2006). 

Using magnetic field data recorded simultaneously by 
the four Cluster spacecraft and assuming that turbu- 
lence contains many structures on scales to be mea- 
sured a nd the time series are at least weakly sta- 
tionary (jPincon fc Lefeuvre 1119911 ). the k-filtering tech- 
nique assumes plane wave geometry and has been ap- 
plied to the magnet osphere and magnetic reconnection 
(ISahraoui et al. ||200 4: Gris on et al~ll2005l: iNarita et al. I 
120051: lEastwood et al. 1 120091: iHuang et al 1 12010^. It is 
well-known that the k-filtering method is subj ect to a 
spatial aliasing effect (jPingon fc Lefeuvre 1 119911 ) . Great 
care must be taken to eliminate or minimize the spatial 
aliasing. This can be realized by setting the maximum 
wave number and spacecraft frequency to be analyzed 
properly (Sahraoui et al. 2010b). Its application to the 
solar wind turbulence is limited and results are inconclu- 
sive: Sahraoui et al. (2010a) found that KAWs populate 
in the solar wind turbulence ion scales while Narita et al. 
(2011) concluded that linear Vlasov theory is insufficient 
to describe the plasma turbulence and turbulent cascade 
is at work. It should be noted that the data studied 
in ISahraoui et al. I (l2010aD were taken during a coronal 
mass ejection ( Jian et al. II2006I ). Narita et al. (2011) 
used data when the tetrahedral configuration of the Clus- 
ter spacecraft was not optimal: the planarity P and elon- 
gation E, which describe the degree that the four Cluster 
spacecraft are close to perfect tetrahedron (Robert et al. 
1998), were such that P > 0.3 and E ^ 0.1, undesirable 
to apply the k-filtering (Sahraoui et al. 2010b) in such 
geometries. 
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In this paper, we present a new study of Cluster data to 
study solar wind plasma turbulence at ion kinetic scales 
by combining the k-filtering technique and wave polariza- 
tion analysis. Although unable to determine wave prop- 
agation direction, polarization analysis supports the in- 
terpretation of KAWs in the turbulence dissipation range 
when interplanetary magnetic field is in the direction 
nearly perpendicular to the solar wind (He et al. 2012). 
We present data analysis in section 3 and discussions on 
the interpretation of the data in section 4. We summarize 
our findings and conclude the paper in section 5. 

2. DATA 

Table 1 summarizes key parameters of four periods 
(PI, P2, P3 and P4) on 31/01/2004 and 29/02/2004 
when the Cluster spacecraft were in the ambient fast so- 
lar wind. The mean parameters of the periods are: Bq 
the strength of the averaged magnetic field, n the to- 
tal ion density, /3 the ion plasma beta (ratio between ion 
parallel thermal pressure and magnetic field pressure) , Vf 
the solar wind speed, fd — e-Bo/(27rmp) proton gyrofre- 
quency, va the Alfven speed, E elongation, P planarity, 
0VB the angle between the solar wind and Bq, T±/T\^ hot 
ion temperature anisotropy, Tg/Ti the ratio of electron to 
ion temperature, pi ion thermal gyro-radius, di = va/^p 
ion inertial length, and Ua/up the abundance of alpha 
particles (fully ionized helium). The rip and are the 
densities of protons and alpha particles. During the cho- 
sen intervals, the magnitude of the magnetic field was 
quite stable and there were no obvious discontinuities 
(see the raw magnetic field data from C4 in P3 in Fig. 
la). Both planarity P and elongation E are smaller than 
0.1 during the periods. 

The magnetic fie ld data were from th e Fluxgate Mag- 
netometer (FGM) (jBalogh et al. 11200 ID . FGM measures 
components of the magnetic field in the GSE (geocentric 
solar ecliptic) coordinate system. In the coordinate sys- 
tem, positive X points from the Earth to the Sun, and 
positive z points to the ecliptic north pole. We use full 
resolution magnetic field data (22 samples/sec). The av- 
erage distance d between the spacecraft was d ^ 200km 
in the four periods. The magnetic field was primarily ori- 
ented in the direction perpendicular to the solar wind di- 
rection so direct magnetic connection with the bow shock 
does not exist. On 31 Jan., the ion pl asma data from the 
Hot Ion Analyzer (HIA) instrument (jReme et al. Il2001| ) 
(with spin resolution) are available from CI spacecraft, 
and on 29 Feb. they are available from both CI and 
C3 (the difference between them is very small). Elec- 
tron temperature data are obtained from the Plasma 
Electron and Current Experiment (PEACE) instrument 
(Johnstone, 1997) onboard the C4. 

3. DATA ANALYSIS 

3.L K-filtering 

Fig. lb shows the Fourier power spectra of the three 
magnetic field components of data from C4 during the 
periods P3 and P4. The spectra are typical of the tur- 
bulent magnetic field fiuctuations in the solar wind. At 
relatively low frequencies (0. 007-0. 4Hz) the fiuctuations 

have an fsc^'^ Kolgomorov power law. At a breakpoint 
fsc ~ 0.4 — 0.5Hz, the spectra steepen with a spectral 
index of about -3.5. The spectra become flattened again 



at the second breakpoin t roughly at 2.4Hz d ue to FGM 
reaching the noise floor (jBalogh et al. |[200l . 

The k-filtering method is a measurement technique de- 
signed for multipoint measurements which does not re- 
quire Taylor's frozen-in fiow hypothesis (Taylor 1938): 
using plane wave assumption, it estimates the spec- 
tral energy density P(w, k) in Fourier space (angular 
frequency lo and wave vector k domains) by combin- 
ing several time series recorded simultaneously at dif- 
ferent locations in space. The k-filt ering method uses 
a fi lter bank approach (jPincon fc M otschniann 199^ 
iTjulin et al] 120051 ) by adopting the random phase ap- 
proximation. The filter is dependent on w and k, and 
is designed in such a way that it absorbs all wave field 
energy except those plane waves with w and k. 

Similar to temporal Fourier analysis, if the spacecraft 
distance is d, the maximum w ave number the space- 
craft can measure is fcma,^ = i i Id (jPincon fc MotschmannI 
IT998I: iSahraoui et al. ll20T0bl) . Due to the use of Fourier 
analysis, spatial aliasing will occur when the spacecraft 
configuration does not distinguish two plane waves dif- 
fering only in wave vectors by Ak: 

Ak • r^j = 27rn,y , l<i<i<N (1) 

where Vij = Vi — Tj is the separation vector between two 
spacecraft (i and j), riij is an integer, and N is the num- 
ber of spacecraft. For Cluster mission (A^=4), the solu- 
tion to the above equation, can be written as (Neubauer 
& Glassmeier 1990, Tjulin et al. 2005): 

Ak ni4Aki + 7124 Ak2 + ?^34Ak3, (2) 

where 

Aki = 27rr24 x v^^a/V, 

Ak2 = 27rr34 x vu/V, (3) 

Aka = 27rri4 x r24/V^, 

and 

y = ri4 • (r24 X r34) 

In our analysis, only wave energy peaks in the k space 
centered at k = are counted by assuming they are due 
to waves physically present in the solar wind (not due to 
aliasing). This k space is given by 

k = ,uiAk+l+^9Ak2-l-/i3Ak3, where -1/2 < //i 2 3 < 1/2. 

' ^ (4) 

A wave with a wave vector in this region will not pro- 
duce any aliased energy peak in the region. Outside 
this region, wave energy peaks will be dropped from our 
analysis. Obviously, a wave with a wave vector outside 
this region may also produce aliased energy peaks inside 
the region. However, this issue is not expected to influ- 
ence our analysis in a significant way due to two rea- 
sons. Firstly, we may generally assume that turbulence 
at smaller wave numbers contains more power than at 
larger wave numbers. Hence, the aliased energy peaks 
produced by larger wave numbers may be too weak to 
be noticeable when the wave energy peaks of small wave 
numbers are present. Secondly, turbulence with larger 
wave numbers may also have higher frequencies (at least 
for normal plasma modes). As a result, the power of 
waves with larger wave numbers may be filtered when 
we analyze the power of waves at low frequencies. We 
notice that when the four Cluster spacecraft form a regu- 
lar tetrahedron configuration, the magnitude of the three 
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wave vectors in Eq. (3) is greater than 2kmax- There- 
fore aliased energy peaks have substantial difference in 
their wave numbers. This fact is strongly in favor of our 
first argument since it is generally known that the so- 
lar wind turbulence power rapidly drops with increasing 
wave numbers. 

The two vertical dashed lines in Fig. 1 represent the 
minimum and maximum frequencies fmin = 0.07Hz and 
fmax=^-^i^z between which the k-filtering technique is 
applied in this paper. To avoid or minimize the spa- 
tial aliasing, a maximum spacecraft frequency fmax has 
to be set corresponding to th e maximum wave number 
kmax (jSahraoui et al. |[2010b[ ). This is also necessary to 
avoid a frequency aliasing effect (Narita et al. 2010). 
Note, it is fortunate that generally high spacecraft fre- 
quency corresponds to large k. In the solar wind rest 
frame, the maximum frequency is kmaxVph / i'^'^) , where 
Vph = Max(v A, Cs) and Cs is the ion sound speed. The 
choice of Vph is equivalent to the assumption that there 
are no whistler waves at scales near kmax in the solar 
wind. This will be verified later on for the solar wind 
data we analyzed. If whistler waves do exist, the choice 
of Vph and the maximum frequency must be dealt with 
accordingly. Since the solar wind is supersonic and su- 
per Alfvenic, the maximum frequency fmax may be set at 
kmax{Vf—Vph)/{2TT) = 1 .32IIz duc to Dopplcr effcct (here 
Vf is the solar wind speed). Note, it is likely that the 
wave vectors will deviate from the solar wind direction 
at an angle 9kVf In such a case, a spacecraft frequency 
higher than fmax = kmax{Vf - Vph) / cos 9 kVj wiU 
correspond to a wave number larger than kmax- In the 
periods we studied, it is found that the vectors can devi- 
ate up to dkVf — 30° from the solar wind direction at the 
highest frequencies (and wave numbers). In this work we 
set fmax — l.lHz. A key reason of choosing this max- 
imum frequency is that above this frequency the noise 
level of FGM is generally believed to be high (according 
to FGM PI Elizabeth Lucek, private communication). In 
fact, within 1.1 < fsc < 1.32IIz, we are still able to use 
the k-filtering if the FGM noise is low. We chose l.lHz 
as the upper limit to avoid producing unphysical results 
due to the FGM noises. 

The fmin value is fixed by choosing kmin = kmax/ "25 
and fmin ~ kminVf / i'^T^) SO that the wave vectors are 
computed with relative ac curacy better than 10% (1% 
at the highest frequency) (jSahraoui et al. ir2010b[ ). Ob- 
viously, the minimum frequency is also limited by the 
number of sampling points available in a dataset. It is 
important to point out a limit of k-filtering technique in 
determining the solar wind turbulence power. For wave 
vectors almost perpendicular to the solar wind direction, 
the wave number component parallel to the solar wind 



flow is fc|| 0. In this case, if the wave has a small fre- 
quency LOpias in the solar wind frame, the Doppler shifted 
frequency fsc = {i-^pias + ^fk\\)/{'2n) may be lower than 
fmin- Such a wavc will not be resolved by k-filtering. 
However, as long as the wave number is in the wave num- 
ber space described by Eq. (4) and the Doppler shifted 
frequency fsc is greater than fmim the wave will be re- 
solved by the k-filtering. 

By scanning the k space, the k-filtering technique is 
used to determine the strongest wave power P{lUsc, k) 
and the corresponding wave vector k at each fsc- The 



wave power in the solar wind frame P{LOpias, k) is then 
determined using the Doppler shift ujpias — ^sc ~ k • V/, 
and the wave dispersion relation ujpias = Wpias(k) is 
obtained. Four studied intervals PI, P2, P3 and P4 
are shown respectively in black, blue, green and red 
in Fig. 2. CIS onboard moments are used for ion 
parameters. A small correction is made to the solar 
wind speed since a few percent of the ions are minor 
ions (mainly fully ionized helium) and the CIS onboard 
moments are calculated by assuming all the detected 
ions ar e protons (HIA/CIS measures the ion energy per 
charge iReme et al. I (j2001[ )). The abundance of helium 
ions can be found from the ion velocity distribution func- 
tion (VDF) measurements by assuming that protons and 
helium ions have the same flow speed Vs so two popula- 
tions can be separated (Marsch et al. 1982). Given the 
solar wind speed from the CIS onboard moment data Vf , 
one finds 



Vs 



1 + ng/up 
1 + 2na/np 



Vf 



(5) 



We use Vs (instead of Vf) to compute ojpias- 

As shown in Fig. 2a, wave vectors are mainly in direc- 
tions quasi-perpendicular t o Bp with (Oun) = 81 ~ 90.1°, 
similar to previous work USahraoui et al. ll2010aD . The 
wave vectors and the solar wind flow make moderate an- 
gles, {OwVf) = 13 — 30°, so generally waves are prop- 
agating in directions not far away from the solar wind 
direction. The error bars (Fig. 2a) at low frequencies 
are significantly larger than at higher frequencies, reflect- 
ing larger relati ve uncertainty when de termining smaller 
wave numbers (jSahraoui et al. |[2010b[ ). 

Fig. 2b displays the measured dispersion relation 
(filled dots). The error bars in the figure mainly come 
from an assumed 3.5% uncertainty in the solar wind 
fiow sp eed in addition to the uncertainty in the wave 
vector (jSahraoui et al. Il2010bf ). At the energy channel 
(2359.28eV) of HIA/CIS where the peak particle flux 
of the solar wind is measured during the studied peri- 
ods, the energy resolution is about < 7%. Hence, the 
error from a 3.5% uncertainty in the solar wind speed 
at 650km/s with Alfven speed at va = 120km/s and 



TABLE 1 

Average plasma parameters during chosen intervals that 
k-filtering technique is applied (data are from cis. fgm 

AND PEACE). 





Jan. 31 (PI) 


Jan. 31 (P2) 


Feb. 29 (P3) 


Feb. 29 (P4) 




14:30-14:40 


14:45-14:55 


04:10-04:20 


04:25-04:35 


B (nT) 


8.45 


7.97 


9.56 


9.34 


n(cm ^) 


3.47 


3.25 


2.88 


2.73 


/3 


0.62 


0.72 


0.73 


0.67 


Vf 


613 


609 


646 


657 


fci 


0.129 


0.122 


0.146 


0.142 


VA 


99.1 


96.2 


123.1 


123.4 


E 


0.05 


0.04 


0.01 


0.02 


p 


0.07 


0.06 


0.03 


0.01 


dvB 


75.1° 


66.6° 


78.6° 


84.1° 




1.41 


1.28 


1.26 


1.46 


Te/T, 


N/A 


N/A 


0.37 


0.39 


pi (km) 


115 


121 


129 


137 


di (km) 


122 


126 


134 


138 




1.4% 


1.3% 


0.38% 


0.2% 
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ikVf 



20° and kvA/^p = 1.4 can be estimated as 



Aw. 



'plas 



kvA ombVf 

X ^ COS 20 



0.25. 



VA 



However error bars in the dispersion plot of Sahraoui et 
aL (2010) are puzzlingly small. Plotted in Fig. 2b are the 
dispersion relations of waves propagating in some mea- 
sured propagation angles, 0kB. They include fast and 
Bernstein waves propagating at 89.5° (red solid and dot- 
dashed), and KAWs propagating at 80°, 85° and 89°. 
At very high k at which our data are unable to cope, the 
branches of KAWs are highly dispersive and are named 
"oblique whistler" waves (Sahraoui et al. 2012). In com- 
puting the dispersion, the abundance of alpha particles 
is 2% and the alphas and protons have the same thermal 
speed. 

Note that ujpias/^p can be negative. It is found that 
95 data points of ujpias are positive while 36 are nega- 
tive. Most of the negative frequencies are small and 27 
of the negative data points are within the uncertainties 
of small positive frequencies. The largest uncertainties 
of ujpias/^p in Fig. 2b are ±0.25 at large kvA/^p ~ 1.4. 
The uncertainties are about ±0.04 when kvA/^p = 0.2. 
We found that 9 of the 36 negative frequencies may have 
to be interpreted as waves propagating in the sunward 
direction in the solar wind frame. Statistically they are 
less important and we will defer their investigation to a 
future study. 

In Fig. 2d, power spectral density of magnetic field 
fluctuations as a function of ky at spacecraft-frame fre- 
quency 0.51Hz in P3, a well-behaved peak, is shown. In 
the X and z directions, the peaks are narrower than in 
the y direction. Therefore, the dominant k is well defined 
in the data. Fig. 2e displays the measured magnetic field 
k± spectra of the four intervals. Few data points exist 
at small k. Data from the four intervals are combined. 
The two solid lines show two power laws with spectral 
indices of -5/3 and -3.5. The spectruin roughly reveals 
two power laws ([Sahraoui et al. Il2010al ): a Kolmogorov 

scaling ^ k^^^^ at smaller k± above a breakpoint at 

k±Pi ~ 0.4 — 0.5. The spectrum steepens to a fcj^'^ 
scaling in an ion dissipation range k±pi e [0.5 — 1.5]. 

3.2. Polarization analysis in the plane perpendicular to 

k 

Once k — (k^^ ky,kz) is found, a primed Cartesian co- 
ordinate system is constructed to study wave polariza- 
tion. The direction of k is along the z'— axis with a unit 
vector e^/ — k/fc. The unit vector along the primed 
x'— axis is e^', and e^i = e^' x eyi. Let 

= [-kyA/k^,A,0] (6) 

describe the three components of e^' in the GSE coordi- 
nates, where A = ^ . Then the three components 

of magnetic field fiuctuations in the GSE coordinates are 
projected on the primed coordinates. A Morlet wavelet 
transform, a natural bandpass filter, is used (He et al. 
2012) and the time series reconstructed at a frequency 
fsc as (Torrence & Compo 1998): 



6Bf 



St^/^ Re{B{fsc)) 



Parameters used for reconstruction are Cg = 0.776, 
tpQ = 7r~4, these are empirically derived for the Morlet 
wavelet. Here, s/^^ (the order of time scale at fsc) is used 
to convert the wavelet transform B to an energy density 
(Torrence & Compo 1998). At each frequency fsc, the 
reconstructed time series contain wave power within a 
frequency window which is about 8.3%/sc centered at 

f sc 

The top panels of Fig. 3 show SB^' — SByi hodograph 
at four frequencies 0.96Hz, 0.74Hz, 0.52Hz and 0.20Hz 
(from left to right) for P3 and P4 using data from C4. 
Results from other three spacecraft are essentially the 
same. Each column corresponds to results of one fre- 
quency. The magnitude of kvA/^p determined by k- 
filtering is 1.35, 0.99, 0.75 and 0.39 for fsc =0.96, 0.74, 
0.52 and 0.2Hz, respectively. The bottom panels show 
the dtp/dt as a function of time, here ip is the angle that 
a magnetic field vector makes with the 5Bx' axis such 
that 



arctan 



5By,{t) 



5B^'{t) 



(8) 



(7) 



Positive (negative) sign indicates that the polarization 
of the wave is right- (left-) handed. It is clear from 
Fig. 3 that the polarization of magnetic fiuctuations in 
the plane perpendicular to the wave vector k is domi- 
nantly right-handed. A closer look of the bottom pan- 
els finds that the dominance of right-handed polariza- 
tion is more pronounced at 0.96Hz, 0.74Hz and 0.52Hz 
than at 0.2Hz. From Fig. 2e, we know that at 0.96Hz, 
0.74Hz and 0.52Hz the wave vectors {kvA/^p is 1.35, 
0.98, and 0.75) are in the dissipation range of the mag- 
netic field power spectrum and at 0.2Hz the wave vector 
{kvA/^p — 0.39) is at (or near) the spectral break point 
where the spectrum switches from inertial range to the 
dissipation range. This may suggest that the turbulence 
has experienced some subtle change in the dissipation 
range where the ion kinetic effect starts to kick in. For 
PI and P2, the polarization is also predominantly right- 
handed at all frequencies studied by the k-filtering tech- 
nique. The Bq in Fig. 3 is the projection of the average 
magnetic field in the plane for the whole period P3 (Figs. 
3c and 3g) or P4 (Figs. 3d and 3h). 

4. INTERPRETATIONS AND DISCUSSIONS 

It is clear from Fig. 2b that the measured disper- 
sion relation cannot be explained by fast or Bernstein 
waves. These measured dispersion relation points are 
quite scattered, and no dispersion relationship of a sin- 
gle plasma wave can be uniquely identified from the mea- 
su red dispersion points , in accordance with the findings 
bv lNar ita et al."] ()2011h . From the k-filtering result, it is 
not clear if we have observed kinetic Alfven waves, con- 
vected coherent structures, or a mixture of them (and 
others) in the solar wind. For instance, while many of 
the data points may be interpreted to be on the dis- 
persion curves of the quasi-perpendicular propagating 
KAWs within the uncertainties, they can equally be said 
to be on the dispersion curve of convected coherent static 
structures within the uncertainties. 

From Fig. 3, except some less frequent anomalies one 
can see that the major axis of magnetic ellipse is dom- 
inantly perpendicular to Bq, and this has been inter- 
preted as evidence of dominant KAWs and not whistler 
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waves (He et al. 2012), although He et al. (2012) have 
to make assumptions on the wave propagation direction. 
During periods P3-P4, the electron temperature is lower 
than the proton temperature (for PI and P2, PEACE 
Te data are not available from the ESA Cluster Active 
Archive) , kinetic slow waves are not expected to exist due 
to strong damping. At 80° < (?kB < 90°, the dispersions 
of fast /whistler waves are all similar to the red-solid line 
in Fig. 2c and too far away from the observed disper- 
sion points. Hence, the wave polarization analysis in the 
plane perpendicular to k may support the interpretation 
that KAWs are an important turbulence component at 
the ion kinetic scale turbulence. 

However, KAWs interpretation has weakness. In the 
studied periods, the polarization is dominantly right- 
handed {{dip/dt), the average dip/dt, is positive) for both 
0kB > 90° and 0kB < 90°. According to Vlasov theory, 
in plasma of one ion species with Maxwellian VDF, the 
magnetic field of KAWs (in the plasma frame) has right- 
handed polarization when 6'kB < 90° and left-handed 
polarization when 6'kB > 90° (here |6'kB - 90°| < 20°). 
The ion plasma betas in this study are smaller than 1 
(used by He et al. 2012) and is only half of T,. We 
find that the change of ion beta (0.6 - 1) and Te/Tp (1 
- 0.5) does not change the polarization of these waves. 
One possibility of the observed {dip/dt) at 6'kB > 90° is 
due to the large uncertainty of 0kB from k-filtering, the 
0kB of the observed left-handed waves is actually smaller 
than 90°. Another weakness of KAWs is that the wave 
power along k is found at least as strong as those in the 
direction parallel to Bq. This is shown in Fig. 4: at 
two frequencies fsc =0.74 and 0.2 Hz, the reconstructed 
fluctuated magnetic field 5Bk and 5B\\ along the direc- 
tion of k and Bq are shown as black and green lines 
within interval P3 (the data are from spacecraft C4). At 
fsc — 0.74Hz, the fluctuated magnetic field along the 
wave vector 5Bk is slightly stronger than the fluctuated 
magnetic field along the background magnetic field. At 
fsc = 0.2Hz, the fiuctuated magnetic field along the wave 
vector 5Bk is often twice as strong as 5B\\ . However, a 
kinetic Alfven wave propagating along a wave vector k 
is expected to generate no fluctuated magnetic field this 
direction. 

Since the measured dispersion points are quite scat- 
tered, they may be seen as no clear dispersion (Narita 
et al. 2011), but the superposition of different things 
such as waves and turbulent structures. An alternative 
interp retation of the d ata is that static small scale cur- 
rents ()Perri et al.ll2012[ ) and 2D nonlinear coherent struc- 
tures (such as solitary monopolar and dipolar Alfven vor- 
tex filaments) with k±^ ^ fc|| populate at ion kinetic 
scales. Monopolar Alfven vortices are static structures 
and dipole Alfven vortices move with an arbitrary speed 
in the plasma frame mainly in the direction perpendicu- 
lar to Bq. The magnetic field fluctuations mainly occur 
in the direction perpendicular to Bq, which is the case 
shown in Fig. 3. Indeed, dispersion relations of these 
static currents and structures are flat in the solar wind 
frame (Fig. 2c). 

To discuss the idea further, we conduct another polar- 
ization analysis in the plane perpendicular to Bq and use 
Bq to replace k in Eq.(6) to construct new Cartesian co- 
ordinates. The results for two frequencies fsc = 0.74Hz 
and 0.2Hz are shown in Fig. 5. On average, the polariza- 



tion of fluctuations {{dip/dt}) can be either positive (Fig. 
5b) or negative (Fig. 5f ) . (The randomness of polariza- 
tion in the plane perpendicular to Bq is fine for KAWs 
since such a wave is supposed to be linearly polarized 
in the plane and the presence of many such waves can 
generate random overall polarization). At each fsc, the 
preference of polarization at one sense is weak {\{dip/dt)\ 
is smaller compared to those in the plane perpendicular 
to k). Magnetic fiuctuations in Figs. 5c and 5g consist of 
wave packets. The hodograms (Figs. 5d and 5h) of the 
perpendicular field 6'B±_ show that Cluster went through 
regions of shear in the magnetic field (labeled as 'C') 
and rotations (for instance at 173.3s, 174.6s and 175.8s 
in Fig. 5d). Such coherent rotations are signatures of 
coherent Alfven vortices (Chemyrev et al. 1988, Volw- 
erk et al. 1996). The rotational sense changes frequently 
(blue and red denote opposite rotations). When the po- 
larization changes (color changes between red and blue) 
in Figs. 5d and 5h, the SB± does not experience any ap- 
preciable change in either magnitude or direction: these 
polarization changes do not correspond to discontinuities 
or currents sheets. 

Alfven vortices (drift Alfven vortices) are 2D tubular 
structures and exist in homogeneous (inhomogeneous) 
plasmas. The observed polarization depends on the tra- 
jectory of satellites across the monopolar or dipolar vor- 
tex: it can be elliptical (linear or circular) , right- or left- 
handed as a function of the trajectory. In homogeneous 
plasmas, there is no limit to the dimension (radius) of 
Alfven vortices. In inhomogeneous plasmas, the theory of 
Alfven vortices valid for scale sizes of ion inertial length 
and ion Larmor radius can be found in Chmyrev et al. 
(1988) and Onishchcnko et al. (2008), respectively. The 
dimension of measured Alfven vortices tends to be the or- 
der of the ion gyroradius (Sundkvist et al. 2005). Such 
drift Alfven vortices are generated naturally in plasmas 
with strong gradients when the drift velocity of parti- 
cles Vd — — Vp X 'B/neB'^ is comparable to their ther- 
mal velocity (Petviashvili & Pokhotelov 1992), or equiv- 
alently the density scale size matches the ion Larmor 
radius (Sundkvist & Bale 2008). 

In the solar wind at lAU the ion drift velocity is small 
due to weak inhomogeneity and Alfven vortices may be 
used to describe these rotational structures. In Fig. 5c 
(0.74Hz), a wave packet typically lasts At =6-8s. The di- 
mension of such wave packet in the direction perpendic- 
ular to Bo is VfAt/ sin{evB) = 3400 5300km, suggest- 
ing that the radius of such structures is a = 13 20.5pi. 
Similar structures with discontinuities have been studied 
in the context of the solar wind (Verkhoglyadova et al. 
2003), and have been found in the Earth's magnetosheath 
(Alexandrova et al. 2006) and Saturn's magnetosheath 
(Alexandrova & Saur 2008). The wave number deter- 
mined by the k-filtering is kvA/^p = 0.98 at 0.74Hz, 
corresponding to a scale of 1160km, in accordance with 
an Alfven vortex with a — 20.5pi (13.5pi) if the vor- 
tex boundary corresponding to the third (second) zero of 
Bessel function of the first kind. Similarly at 0.2Hz (Fig. 
5g), a wave packet typically lasts 20-22s. The dimension 
of such a wave packet in the direction perpendicular to 
Bo is VfAt/sm{evB) = 13200 - 14500km. The radius 
of such a structure is a = 96 ^ 106pi, considerably larger 
than that at higher frequency 0.74Hz as we would expect. 

The theory of solitary Alfven vortex is based on single- 
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fluid MHD and assiimos incomprcssibility (Pctviashvili 
& Pokliotelov 1992). At the scales wo studied, it is 
clear from Fig. 3 that 5B\^ <C 6B± so the turbulence 
incompressibility is approximately met. The limitation 
of solitary Alfven vortex interpretation is the difficulty 
to explain the polarization in the plane perpendicular to 
k in Fig. 3 with solitary Alfven vortices. This is be- 
cause that the theory of solitary Alfven vortex assumes 
(5S|| ?a 0. A nonzero 6B^^ is necessary to explain the 
dominantly right-handed polarization in the plane per- 
pendicular to k if the wave vector is perfectly perpendic- 
ular to Bq. The k- filtering analysis finds that the wave 
vector mainly points to directions nearly perpendicular 
to Bq. One would expect that the wave polarization in 
the plane perpendicular to k does not have a preference 
in either left-handed or right-handed sense when a space- 
craft passes through many of such structures. A theory 
of solitary Alfven vortex including small compressibility 
is needed. 

5. SUMMARY AND CONCLUSION 

In summary, the application of the k-filtering tech- 
nique and wave polarization analysis to turbulence at 

the proton gyroscalcs in the fast solar wind found the 
following: Turbulence at these scales slowly (compared 
to the Alfven speed) propagates in the directions nearly 
perpendicular to Bq. The fluctuated magnetic field in 
the frequency range 0.07-l.lHz shows higher SB± than 
and has dominantly right-handed polarization in the 
plane perpendicular to the wave propagation direction 
at both ^kB < 90° and ^kB > 90°. The polarization of 
the fluctuations is elliptical with a regular change of po- 



larization from right to left-handed. Wave polarization 
is quite random in the plane perpendicular to the back- 
ground magnetic field and is consistent with the interpre- 
tation of Alfven vortices. The wave polarization in the 
plane perpendicular to wave vector k is more consistent 
with linear kinetic Alfven waves than Alfven vortices. It 
is found that no dispersion relation of a single plasma 
wave mode can be uniquely identified from the mcasiired 
wave/turbulence dispersion plots. 

We have discussed the pros and cons of KAWs and co- 
herent structures in the interpretation of the solar wind 
turbulence at ion kinetic scales. A plausible scenario is 
that at such scales KAWs and coherent structures co- 
exist in the fast solar wind described in this study. It is 
noted that further validation of the k-flltering technique 
may be needed when the analyzed signal contains a mix- 
ture of coherent structures and plane waves with random 
phases, not just plane waves with random phases alone. 
On the other hand, one may see certain similarity be- 
tween a series of intermittent coherent structures of sim- 
ilar sizes passing a spacecraft and a plane wave with a 
random phase passing a spacecraft. We plan to publish 
such a validation elsewhere. 
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Fig. 1.— Spectra of typical magnetic field components B^, By and B^ measured by FGM from 04:10UT to 04:35UT on 29/02/2004. 
The vertical dashed lines denote the frequency range that k-filtering technique is applied. The spectral flattening above 2.4Hz is due to 
the FGM reaching the noise floor. 




Fig. 2. — (a) Angles 6kB (squares) and 9kVj: (solid dots) with related uncertainties computed by using the k-filtering technique during 
four time intervals, (b) Measured wave dispersion (filled dots), with estimated error bars, (c) The dispersion relation curves, computed 
from linear Vlasov theory, represent waves propagating at several observed angles SkB- black lines are kinetic Alfven waves propagating 
at 80° (dashed), 85° (dotted) and 89° (solid); the remaining red (solid and dotted) curves represent fast and Bernstein waves propagating 
at 89.5°. The proton angular gyrofrequency is Qp = 2iTfci. The dot-dashed line represent static structures. (d)Power spectral density of 
magnetic field fluctuations as a function of ky at spacecraft-frame frequency 0.51Hz in P3. (e) The magnetic field k± spectra of all the four 
measured time intervals. 



Kinetic turbulence in the fast solar wind 




Fig. 3. — Polarization analysis in the plane perpendicular to k. Top panels: SB^i — SByi hodograph at frequencies (a) fsc = 0.96Hz, 
(b)0.74Hz in P3, and (c) 0.52Hz, (d) 0.2Hz in P4. The bottom panels display the corresponding polarizations at these frequencies. The 
wave propagation angle S^g is a) 92.6°, b) 95.2°, c) 84.6°, and d) 87.2°. The mean values of [dip/dt) are all positive. The wave frequencies 
in the solar wind frame determined by the k-filtering are all positive. 
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Fig. 4. — Reconstructed time series of fluctuated magnetic field along wave vectors found by the k-filtering technique and along the 
background magnetic field Bq for the interval P3 at two spacecraft frequencies: (a) fgc = 0.74Hz, and (b) fgc = 0.2Hz. The data are taken 
from spacecraft C4. 
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Fig. 5. — Polarization analysis in the plane perpendicular to Bq at /sc = 0.72Hz (a,b,c,d) and 0.21Hz (e,f,g,h) for interval P3. (a) and 
(e); bB^^i —SB^yi hodograph, (b) and (f): dLp/dt, (c) and (g): representative waveforms, magnetic field hodograms for the region between 
the two vertical lines are shown in (d) and (h). Blue and red denote right- and left-handed polarization. At the two frequencies, the wave 
vectors determined by k-filtering are fcii^/tJo =0.981 and 0.354. The data are taken from spacecraft C4. 



